Priority diffusion model in lattices and complex networks 
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We introduce a model for diffusion of two classes of particles {A and B) with priority: where 
both species are present in the same site the motion of A's takes precedence over that of _B's. This 
describes realistic situations in wireless and communication networks. In regular lattices the diffusion 
of the two species is normal but the B particles are significantly slower, due to the presence of the 
A particles. From the fraction of sites where the B particles can move freely, which we compute 
analytically, we derive the diffusion coefficients of the two species. In heterogeneous networks the 
fraction of sites where B is free decreases exponentially with the degree of the sites. This, coupled 
with accumulation of particles in high-degree nodes leads to trapping of the low priority particles 
in scale-free networks. 



Diffusion, or the random motion of particles is a most 
basic mechanism underlying numerous phenomena in na- 
ture and technology. While diffusion in periodic regular 
lattices is rather simple, it is significantly richer in disor- 
dered media and complex networks, or when the particles 
interact with one another. In most studies the particles 
interact through the excluded volume effect, combined 
with a potential (as in 'lattice gas' models), or through 
some kind of reaction or transformation of the particles 

(e.g. aaii). 

In this paper we introduce a two-species priority dif- 
fusion model (PDM): the two species, A and B, diffuse 
independently, but where both species coexist only the 
high priority particles, A, are allowed to move. This 
problem has several important applications. A frequent 
case in communication networks is that data packets tra- 
verse the networks in a random fashion (e.g. in wireless 
sensor networks_[5i, iQ] , ad- hoc networks [7,, ^8,] and peer-to- 
peer networks (9|). Routers in communication networks 
handle both high and low priority information packets, 
such as, for example, in typical multimedia applications. 
The low priority packets are sent out only after all high 
priority packets have been sent [loL [ll|, just as in our 
model. 

We solve the priority diffusion model (PDM) analyt- 
ically for lattices and networks. In lattices and regular 
graphs both species diffuse in the usual fashion, but the 
low priority B's diffuse slower than the j4's. In heteroge- 
neous scale-free networks the i3's get mired in the high 
degree nodes, effectively arresting their progress. We 
confirm these conclusions through large-scale computer 
simulations. 

The A and B particles, when selected for motion, hop 
to one of the nearest neighbor sites, with equal proba- 
bility. We have investigated two selection protocols. In 
the site protocol a site is selected at random: if it con- 
tains both A and B particles, a high-priority A particle 
moves out of the site. A particle of type B moves only 
if there are no A's on the site. If the site is empty, a 



new choice is made. In the particle protocol a particle 
is randomly selected: if the particle is an A it then hops 
out. A selected B hops only if there are no A particles on 
its site 112]. The site protocol describes the case where a 
selected router sends out packets of information, whereas 
in the particle protocol agents move independently, de- 
scribing perhaps the commuting of various individuals. 
Note that these protocols belong to the general frame- 
work of zero-range process with two species of particles 
(see e.g. [ill, [13] for factorized steady-state solutions). 

We note that the A particles essentially move freely (in 
both protocols), regardless of the i?'s. The i?'s, on the 
other hand, can move only in those sites that are empty 
of ^'s. We begin by considering the number of such sites. 
We later relate this property to diffusion coefficients of 
the particles under the priority constraints. 

We look first at lattices, or regular graphs, where each 
site has exactly z nearest neighbors. The number of sites 
— > oo and for now we focus on a single species, denot- 
ing its particle density by p. Let fj be the average (equi- 
librium) fraction of sites that contain j particles. Con- 
sider a Markov chain process whose states {0, 1, 2, ...} are 
the number of particles in a given site. The {/j}j=o,i,2,... 
are the stationary probabilities of the chain. 

For the site protocol, the transition probabilities are: 

= ^ ; ^w+i = • (1) 

Pj^j = 1 — Pj.j_i — Pjj+i and all other transitions cannot 
occur. Indeed, for a site to lose a particle it needs to be 
selected, with probability To gain a particle, a non- 
empty, one of its z neighbors must be chosen — with 
probability (1 — /o) — and this neighbor must send the 
particle into the original site, with probability K Note 
that the final result is independent of the coordination 
number z. The stationary state satisfies 
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with the boundary condition /i = (1— /o)/o- This has the 
solution fj = /o(l — /o)''- Imposing particle conservation 
J2T=o i/i = ^'^ finally obtain: 
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This agrees nicely with our simulation results (Fig. ©a), 
inset). 

For the particle protocol the transition probabilities 
are: 
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Pj j = 1 — Pj.j-i — Pj.j+i and all other transitions are 
excluded. Indeed, for a site to lose a particle one of its 
j particles (out of the total Np) needs to be selected. 
To gain a particle, one of the zp particles that reside, 
on average, in the neighboring sites has to be chosen, 
and then hop to the original site (with probability i). 
Once again, the result is independent of z. This time the 
boundary condition is /i = pfo, leading to fj = /o jr- 
Imposing the normalization condition j^o J"] ~ "^"^ 
finally finds 
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In other words, the {fj} are Poisson-distributed, with 
average p. 

We now employ these results for the analysis of prior- 
ity diffusion, when both species are involved. In regular 
graphs both species diffuse as in the single-species case, 
but due to the priority constraints the total time available 
is split unevenly between the ^'s and B's. For example, 
in lattices diffusion is normal, (R^) — Dt, as confirmed 
by simulations (Fig. [IJa)), but with a smaller diffusion 
coefficient for the hindered B's. 

Denote by Pa (resp. Pb) the probability that a moving 
particle is an A (B). In the site protocol, a particle will 
surely move if we choose a non-empty site (contains A, 
B or both), which happens with probability 1 — 1/(1 -f 
PA + Pb) (since the particles behave as a single, non- 
interacting species, if one ignores their labelling, and thus 
eq. ([U can be applied with p ~ pa + Pb)- A moves if the 
selected site contains any number of ^'s, which happens 
with probability pa/ {^ + Pa): again, from Therefore, 
^^ = Tf^/(l---^)^or: 
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where ps = Pa + Pb and we have used Pb = 1 — Pa 
for the second relation. For a single particle in a lat- 
tice, (-R^) = t- Here we have (i?^) = PAt = DAt and 
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FIG. 1: (Color online) PDM in lattices (site protocol). Typ- 
ical system size is 10* sites and time is 10"^ steps, (a) Mean 
square displacement of a particle (R^) (averaged over all par- 
ticles and realizations) as a function of time. In one time step 
each particle moves on average once. Black squares are for 
densities pA ~ Pb = 10, green triangles for pA = Pb ~ 1- 
Inset: /o as a function of density. The solid curve represents 
Eq. Q, while symbols denote simulation results, (b) Dif- 
fusion coefficients, D {Pa full symbols, Pb empty symbols), 
as a function of pA- Black squares are for pA = lOps, red 
circles for pA = Pb and green triangles for pA ~ O.lps. The 
continuous (Pa) and dotted (Pb) curves represent Eq. (0. 



thus the diffusion coefficient Da equals Pa (similarly, 
Db — Pb)- Simulations shown in Fig. [T)d confirm our 
predictions (Eq. ([7])). 

For the particle protocol, denote the ratio of free B 
particles (that do not share a site with ^'s) to all B par- 
ticles by r. Any particle will surely move except for the 
case when a non-free B particle was chosen (which hap- 
pens with probability pb/ Ps ■ (1 ~^))- P moves whenever 
a free B is chosen, with probability pb/ps ■ r. Therefore 
we find: 
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and Pa = 1 - Pb. 

Had the density of -B's been independent of the A's 
then r would simply be the fraction of sites empty of A, 
or r = e'P^. However, due to the priority constraints 
B's tend to stick with the A's, so that r < e"''-*. The 
ratio r can be obtained analytically for low densities, if we 
assume that a single site cannot contain more than one A 
or B. We use again a Markov chain formulation, but now 
with just four possible states to each site: {(j), A, B, AB} 
(state A corresponds to a site having one A particle, and 
similarly for the other states). We write the transition 
probabilities as before, to first order in the densities: 
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Unindicated transition probabilities are zero, and the di- 
agonal accounts for normalization P^.x = 1 — '^y^x ^^^v 
The justification is similar to that of Eq. ([5]). For a site 
to lose a particle, this particle needs to be chosen out of a 
total of N ps particles. For a site to gain an A, one of the 
zpA particles that reside, on average, in the neighboring 
sites has to be chosen (out of Nps), and then sent to the 
target site, with probability ^ (and likewise for gaining 
a B). The priority constraint is taken into account by 
forbidding the transition AB A. 

From the stationary probabilities of the chain (Eq. ([2|)) 
we derive r to first order: 
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(p stands for either p a or pb)- To obtain the next order, 
allowed states can have two particles of each type, and we 
take into account that when a i? is chosen it actually hops 
only with probability r (using its first-order expression, 
Eq. (Uni)). We thus find 

r = l-2pA + ^pI + Oip^) . (11) 



Surprisingly, r does not depend on ps, at least to sec- 
ond order. In fact, our simulations suggest that r is in- 
dependent of Pb for all densities (Fig. [IKa)). For large 
Pa we have r — > e"''^, since free B's become extremely 
rare. The diffusion coefficients obtained on substituting 
Eq. ()lip in ([8]) compare quite nicely with simulations 
(Fig.Hb)). 

We now turn to heterogeneous networks, where the 
degree k varies from site to site. We focus on two ex- 
plicit models: Erdos-Renyi (ER) random graphs [ll,[lH, 
where the degrees of the nodes are narrowly (Poisson) 
distributed, and scale-free (SF) networks, recently dis- 
covered to best describe a wide variety of natural and 
man-made systems, and in particular many communica- 
tion networks such as the Internet. In SF nets the degree 
distribution is broad, characterized by a power-law tail 
P{k) - k-y, when usually 2<7<3[l3,[li,IIl,[2i. 

We analyze the particle protocol only — the site proto- 
col yields similar results, as we confirmed through com- 
puter simulations. We wish to find the fraction of empty 
sites of degree k, /q . As before, consider a network 
with only one particle species and define a Markov chain 
on the states {0, 1,2...} for the number of particles in a 
given site of degree k. The stationary probabilities are 
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It can be shown that the chain has the transition 



probabilities: 
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Pjj = 1 — Pj.j~i — Pj,j+i and all other probabilities are 
zero. Intuitively, this is the same as Eq. (O, except that 




FIG. 2: (Color online) PDM in lattices (particle protocol). 
Simulation parameters are same as in Fig. [T] (a) Fraction r 
of free B's as a function of pA for the cases pA = Pb (black 
squares), pA ~ 0.2pB (empty red circles), and pA = 5ps 
(green triangles), compared to the result of Eq. (|lip (solid 
line), r is smaller than e~''^ (dotted red line) since B parti- 
cles tend to stick with the ^'s. Inset: Large pA regime, where 
r tends to e~''^ {pB ~ Pa)- (b) Diffusion coefficients (Pa 
full symbols, Pb empty symbols) as a function of pA for the 
cases Pa = lOps (black squares), pA = Pb (red circles), and 
Pa ~ O.lpfl (green triangles). The continuous (Pa) and dot- 
ted (Ps) curves represent Eqs. ((H)) and Inset shows Pb 
(red circles) for large pA- The B particles diffuse much slower 
here compared to the site protocol (continuous red line) due 
to the exponential decrease in the number of empty sites. 



here a site may gain a particle from any of its k neighbors, 
and the neighbor delivers the particle to the target site 
with probability {(k) is the average degree of the net). 
Solving for the stationary probabilities while keeping in 



mind that ^ /^''^ = 1 one finds 
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Note that for regular graphs, when all sites have the 
same degree, this reduces to Eq. fo = e^P. Also, 
the total number of particles in a site of degree k is 
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oc k, as is well known for ran- 
For the full network, 
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In particular, for ER graphs, 
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The agreement between these predictions and simula- 
tions is shown in Fig.[3]Ja). 

Consider now priority diffusion. Define that in one 
time step each particle has on average one moving at- 
tempt. On average, at every time step a B particle (in a 
node of degree k) has a probability exp{~pAk/ {k)) to be 
able to jump out (Eq. (jl3p V This results in a distribution 
of waiting times (for a B particle) 
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FIG. 3: (Color online) PDM in networks. Typical system size 
for both ER and SF is 10* nodes, Internet is 21955 nodes 
and time is 10* MCS. SF networks were generated using the 
MoUoy-Reed algorithm 23]. (a) Fraction of empty sites in 
ER (empty symbols) and SF (full symbols) networks as a 
function of site degree. The curves follow Eq. (|13p . Inset: 
Fraction of empty sites for ER vs. (fc) for fixed density p = 1. 
The continuous curve represents Eq. (|14|l . (b) Distribution 
of waiting times for B particles for nodes of various degrees 
{k = 2, 4, 7, 10, 15, left to right) in SF networks (7 = 3, m = 
2). Inset: (t) as a function of k. Shown are: SF and ER 
networks with (k) = 4 (black squares and green triangles), 
and (fc) = 3 (red circles and blue triangles), and the Internet 
[23 | (cyan diamonds). Solid curves represent the theoretical 
{t{k)) — exp{pAk/ (k)). Plateau for large k results from finite 
simulation time. 



where r = {t{k)) ~ exp(pAk/ (k)) is simply the inverse of 
the probability for the site to be empty of ^'s (Fig.[3I^b)). 

The exponentially long waiting time (in the degree 
k) means that in heterogeneous networks such as scale- 
free nets — where the degrees may span several orders 
of magnitude — the B particles get mired in the hubs 
(high degree nodes). The problem is exacerbated by the 
fact that the B particles are drawn to the hubs even 
in the absence of A's: the presence of A's only amplifies 
this tendency, because of the positive correlation between 
the concentrations of the two species. Thus, while the 
concentration of the A's is proportional to k, it can be 
shown that the concentration of the B's is proportional 
to kexp{pAk/ (k)). In large scale-free nets the B's col- 
lect at the hubs, which results in extremely long waiting 
times and practical halting of particles diffusion. 

In conclusion, we have introduced and analyzed a 
model of priority diffusion, involving two species, where 
the particles of the preferred species always move ahead 
of the other. In regular graphs and lattices diffusion of 
each of the species is similar to that of a single (non- 
interacting) species on the same substrate, but the over- 
all diffusion rate of the preferred species is significantly 
faster: we have provided exact expressions for the rate 
ratios in this case. In heterogeneous nets, we have shown 
that the slow-species particles tend to be collected in the 
high degree nodes, where their waiting time for clearing 
the site increases exponentially with the degree of the 



site. In scale-free nets where the degrees span several or- 
ders of magnitude the slow species progress is effectively 
arrested. 

Our calculations for heterogeneous networks are mean- 
field in character, in that we ignore possible correlations 
between the degrees of neighboring nodes. This assump- 
tion is valid, though, for ER random graphs and for max- 
imally random networks, such as MoUoy-Reed SF nets 
[^3] . The possible effect of degree-degree correlations re- 
mains a subject for future research. 
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